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Abstract 

Physically speaking, the delta function like beam-beam nonlinear 
forces at interaction points (IPs) act as a sum of delta function nonlin- 
ear multipoles. By applying the general theory established in ref. |IJ, 
in this paper we investigate analytically the beam-beam interaction 
limited dynamic apertures and the corresponding beam lifetimes for 
both the round and the flat beams. Relations between the beam-beam 
limited beam lifetimes and the beam-beam tune shifts are established, 
which show clearly why experimentally one has always a maximum 
beam-beam tune shift, £,y : max, around 0.045 for e + e~ circular colliders, 
and why one can use round beams to double this value approximately. 
Comparisons with some machine parameters are given. Finally, we 
discuss the mechanism of the luminosity reduction due to a definite 
collision crossing angle. 

1 Introduction 

Beam-beam interactions in circular colliders have many influences on the 
performance of the machines, and the most important effect is that beam- 
beam interactions contribute to the limitations on dynamic apertures and 
beam lifetimes. Due to the importance of this subject, enormous efforts 
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have been made to calculate incoherent and coherent beam-beam forces, to 
simulate beam-beam effects, to find the difference between flat and round 
colliding beams, and to establish analytical formulae to estimate the maxi- 
mum beam-beam tune shift P|-||20||. Physically speaking, the delta function 
like beam-beam nonlinear forces at interaction points (IPs) act as a sum of 
delta function nonlinear multipoles. In ref. Jl| we have established a gen- 
eral theory to study analytically in detail the delta function multipoles and 
their combined effects on the dynamic apertures in circular storage rings, and 
in this paper we will apply these general analytical formulae to the case of 
beam-beam interactions and find the corresponding beam dynamic apertures 
and beam lifetimes. We will show quantitatively why there exists a maximum 
beam-beam tune shift, £, y ,max, around 0.045 for flat beams in e + e~ circular 
colliders, and why this number can be almost doubled for round colliding 
beams. Applications to some machine parameters are also given. In this 
paper we will restrict ourselves to the discussion of e + e~ circular colliders 
since the treatment for the hadron colliders will be somewhat different and 
more difficult. Finally, we discuss beam-beam effects with a definite crossing 
angle. 



2 Beam-beam interactions 



For two head-on colliding bunches, the incoherent kick felt by each particle 



can be calculated as [13 



Sy' + iSx' 



N e r e 
7* 



f(x,y,a x ,a y ) 



(1) 



where x' and y 1 are the horizontal and vertical slopes, N e is the particle pop- 
ulation in the bunch, r e is the electron classical radius (2.818 xlO -15 m), a x 
and Uy are the standard deviations of the transverse charge density distri- 
bution of the counter-rotating bunch at IP, 7* is the normalized particle's 
energy, and * denotes the test particle and the bunch to which the test par- 
ticle belongs. When the bunch is Gaussian f(x,y,a x ,a y ) can be expressed 



by Basseti-Erskine formula | 1Q| : 



f(x,y,a x ,a y 
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where w is the complex error function expressed as 



it? (z) = exp(— z 2 )(l — erf (— iz)) (3) 

For the round beam (RB) and the flat beam (FB) cases one has the incoherent 
beam-beam kicks expressed as [|] [p!2[p3| : 

(RB, a x = a y = a) (4) 

exp(u 2 )du (FB, a x >> a y ) (5) 
(FB, a x » (j y ) (6) 




where r = \Jx 2 + y 2 . Now we want to calculate the average kick felt by the 
test particle since the probability to find the transverse displacement of the 
test particle is not constant (in fact, the probability function is the same as 
the charge distribution of the bunch to which the test particle belongs in 
lepton machines due to synchrotron radiations). In the following we assume 
that the transverse sizes for the two colliding bunches at IP are exactly the 
same. For the round beam case after averaging one gets0][I3|: 



^-^?(l-exp (-£!)) (RB, (7) 

Although this expression is the same as that of the coherent beam-beam kick 
for round beams, one should keep in mind that we are not finding coherent 
beam-beam kick originally, and the difference will be obvious when we treat 
the vertical motion in the case of flat beams. For the flat beam case, we will 
treat the horizontal and vertical planes separately. As far as the horizon- 
tal kick is concerned, the horizontal kick depends only on one displacement 
variable just similar to the round beam case, we will use its coherent form 
expressed as follows 



where a x in the incoherent formula in ref. [|T^| has been replaced by T, x = 
\plo x (for two identical Gaussian colliding beams) according to Hirata the- 
orem demonstrated in the appendix A of ref. |14| . As for the vertical kick, 
however, one has to make an average over eq. ^| with the horizontal proba- 



bility distribution function of the test particle, and one gets |T3 



where <> x means the average over the horizontal probability distribution 
function of the test particle, and for two identical colliding Gaussian beams 
<> x = 1/V2- It is obvious that eq. |9] is not the expression for the coherent 
beam-beam kick. The average over eqs. |4] and |6] is only a technical operation 
to simplify (or to make equivalence) a two dimensional problem to a one 
dimensional one. To study both round and flat beam cases, we expand Sf' 
at x = (for round beam we study only vertical plane since the formalism 
in the horizontal plane is the same), 5x' and 5y' expressed in eqs. [7], || and 
H respectively, into Taylor series: 

N e r e .1 1 o 1 ^ 1 7 

5 = — — f v v H v v 

y 7* l 2a 2 ^ 16a 4 y 192a 6 y 3072a 8 y 

+ 7^r-my 9 - n^L ^ n + »Jo n i, y 13 ----) ( RB ) ( 10 ) 



61440a 10 1474560a 12 41287680a 
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o„ = -x H -x -x 

2 7# V a 2 3a* 30a 6 420a 8 

1 1 11 1 13 

+ 7560a 10 V ~ 166320a 12 % + 4324320a 14 % ) (FB) (H) 

r/ N eT e ,2 1 3 1 5 1 - 

V27* 3a x aj; 20a a: a° 1680-^ 

+ l —J TTV 11 + " T^V 13 -■■■) (FB) (12) 

1728a ;E a9 y 21120a z a 11 y 299520a x a 13 y ' v ; v ' 

The differential equations of the motion of the test particle in the transverse 
planes can be expressed as: 
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1 00 

+ 299520^^- '■') ^ (FB) (15) 

where K x (s) and K y (s) describe the linear focusing of the lattice in the hor- 
izontal and vertical planes. The corresponding Hamiltonians are expressed 
as: 

2 2 y 7* l 4a 2i/ 64a 4i/ 1152a 6y 

1 oo 

' 2/ 8 + ---) E ^( s — kL) (RB) (16) 
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1344a ;E a„ 
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where p x = dx/ds and p y = dy/ds. 
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3 Review of the general analytical formulae 
for dynamic apertures 

In ref. |l| we have studied analytically the one dimensional (y = 0) dynamic 
aperture of a storage ring described by the following Hamiltonian: 

rr P 2 K (s) 2 1 d 2 B z . ~ 

H = — + -^x 2 + — : 3 L Y 6(s - kL) 

2 2 3\Bp dx 2 fc f^ 1 } 

1 o3 TO oo 

4!5p dx 3 k ^ 

where 

fi 2 = B (l + xb x + x 2 6 2 + x% + x% + ■■■ + x m - x b m _ x + ■■■) (20) 
The dynamic aperture corresponding to each multipole is given as: 

l/(m-2) 





2(m-2) 
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\mf3™(s 2m )j 
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A dyna>2m , x (s) = J2f3 x (s) —-j- . (21) 



where S2 m is the location of the 2mth multipole, f3 x {s) is the beta function in 
x plane. Since these results are general, we have tried to avoid assigning the 
freedom of motion, x, a specific name, such as horizontal, or vertical plane. 



4 Beam-beam limited dynamic apertures 

To make use of the general dynamic aperture formulae recalled in section 3, 
one needs only to find the equivalence relations by comparing three Hamilto- 
nians expressed in eqs. [H], [T7|, and [TE| with eq. [H5], and it is found by analogy 
that: 

bm-i T N e r e 



(RB) (22) 
L= ^ ^ e 'L„ m (FB,x) (23) 

(FB, y) (24) 



bm-i L _ N e r, 
b m -\ N e r e 
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12 


14 




16 


192 


3072 


61440 


1474560 


41287680 


C m ,FB,x 


3 


30 


420 


7560 


166320 


4324320 


C mt FB,y 


3 


20 


168 


1728 


21120 


299520 



Table 1: summary of multipole coefficients 
where C m ,RB, Cm,FB,xi an d C m ^FB,y are given in Table 1. Now by inserting 



eqs. |22|-|24| into eq. |21] one can calculate dynamic apertures of different 
multipoles due to nonlinear beam-beam forces. For example, one can get the 
dynamic apertures due to the beam-beam octupole nonlinear force: 



>-dyna,8,y 




(RB) (25) 



(26) 



Pv(sip) \ N e r e 



(FB) (27) 



where sjp is the IP position. Given the dynamic aperture of the ring without 
the beam-beam effect as A x y , the total dynamic aperture including the beam- 
beam effect can be estimated usually as: 

A to tai,x,y{s) = i 1 1 === (28) 
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In the following we will consider the case of A to tai,x,y(s) ~ Abb )X ,y(s)- If we 
measure the beam-beam interaction limited dynamic apertures by the beam 
sizes (the normalized beam-beam limited dynamic aperture), one gets: 



71 
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y,8 
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N e r e (3 y (s IP ) 
fry*** 



Adyna,8,y \&) 
<7*M 



\N e r e (3 x (s IP ) 
\N e r e (3 y (s IP ) 



1/2 



1/2 



1/2 



(RB) 
(FB) 
(FB) 



(29) 
(30) 
(31) 



Recalling and using the definitions of the beam-beam tune shifts £ x and £ y 

N e r e P XjIP 



in eqs. [32| and [33 



S.t 



2-Kj*a x (a x + a y ) 
N e r e P y ,ip 



(32) 
(33) 



2^7*0-^(0-;,. + (Jy) 

one can simplify the above defined normalized dynamic apertures. As general 
results one finds: 
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1 

m-2 



1 

m — 2 



1 

m — 2 



(RB) 



(FB) 



(FB) 



(34) 



(35) 



(36) 



Obviously, the normalized beam-beam effect limited dynamic apertures are 
determined only by the beam-beam tune shifts. The impact of this discovery 
will be more appreciated later. When the higher order multipoles effects 
(2m > 8) can be neglected eqs. |25|, ^ and |27| give very good approximations 
dynamic apertures limited by one beam-beam IP. If there are Njp interaction 
points in a ring the dynamic apertures described in eqs. ^ and |27| will be 
reduced by a factor of ^/Njp (if these Njp interaction points can be regarded 
as independent). 
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5 Beam lifetime due to beam-beam 
interactions 

We take the beam-beam limited dynamic aperture as the rigid mechanical 
boundary, i.e., those particles which walk beyond this virtual boundary will 
be regarded lost instantaneously. Based on this physical point of view we 
can calculate the beam lifetime due to the nonlinear beam-beam effect: 



Ub == — 



a y {s) 



exp 



^■dyna,y (^) 

, a y( s ) 2 , 



(37) 



where r y is the synchrotron radiation damping time in vertical plane. It is 
interesting to note that eq. ^ is similar to but different from the quantum 
lifetime forluma where that a y (s) 2 is used instead of 2a y (s) 2 . The reason 
is that the quantum radiation results in energy fluctuations of an electron and 
2a y (s) 2 corresponds to the average of the square of the oscillation amplitude, 
however, the dynamic apertures calculated above due to nonlinear forces 
are relevant only to the projected motions. When the beam-beam octupole 
nonlinear force dominates the dynamic aperture, by inserting eqs. E9l |30| 



and 31 into eq. 37, or inserting eqs. C34L 35, and pq into eq. 37 one gets 




167*0^ 



N e r e (3 y (s IP ) ) J \N e r e p y (s IP 



r bb,x — 



67, 



N e r e (5 x (s IP 




167*a z 



67, 



or 



gxp I - - 

- x / J \N e r e f3 x (sip) 

t* ( 3V2^a x a y \ 1 / 3y/2^a x a y 
Tbh >y~ 2 \N e r e (3 y {s IP )j eXP \N e r e (3 y (s IP ) 



(RB) 

(FB) 
(FB) 




(RB) 
(FB) 



eXP fej (FB) 



(38) 

(39) 
(40) 

(41) 
(42) 
(43) 
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More generally, one has: 



T, 



T * / 2 2 C m j;j; 



2 / 2 \ 

m—2 / / - m-2 \ m _2 

^ 2 ^m,RB 



bb,2m,y 



\ 



47rVm£* 



r 



r, 



bb,2m,x 



2^a 



2 

' m-2 



m,FB,x 



7r2 v /m£* 



2^G 



m,FB,x 



T, 



bb,2m,y 



2 



2_ / 

m—2 I 



m,FB,y 



Ti2y/mi* 
7T\/2m£* 



2 

m-2 



2 

m—2 



exp I - I I (RB) (44) 

exp | _ n I I (FB) 

(45) 

exp | — ' ' | | (FB) 

(46) 

If we define the lifetime divided by the corresponding damping time as nor- 
malized beam lifetime, one finds that the beam-beam effect limited normalized 
lifetimes depend only on beam-beam tune shifts. Figs. |ll| and |12| show the 
normalized beam lifetime with respect to the beam-beam tune shifts for both 
flat and round beams. 



6 The maximum beam-beam tune shifts for 
flat and round beams 

Now it is high time for us to discuss the maximum beam-beam tune shift 
problem. In literatures the term "maximum beam-beam tune shift" of a 
specific machine is not well defined. One of the reasonable definitions would 
be that the maximum beam-beam tune shift corresponding to a well defined 
minimum beam-beam limited lifetime. In this paper we propose to take this 
well defined minimum beam-beam limited lifetime as one hour (the idea is 
to reduce eq. [281 to A tota i(s) ~ A bb (s), and to have a machine still working !). 



Assuming that for both round and flat beam cases one has the same r y , from 



eqs. 41, 42 and 43 one finds the following relations: 



,rb = V2 1 fb = 1 m FB (47) 

^y,max g >>y,max max \ / 

and 

pFB _ x /2f FB (AS.) 
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It is proved theoretically why round beam scheme can almost double the 
£y,max of flat beam scheme as previously discovered in the numerical simu- 
lations [IB] [17], and why the vertical beam-beam tune shift reaches its limit 



earlier than the horizontal one. Quantitatively, taking r y = 30 ms, one finds 
that Cy,max,FB(nb = 1 hour)= 0.0447, £ x>maXtF B(nb = 1 hour)= 0.0632, and 
i y ,max,RB{nb = 1 hour)= 0.0843. 

Now we investigate how the order of nonlinear resonance affects the max- 



imum beam-beam tune shift. By using eqs. [H], |42|, and [If], and assuming 
that t x = T y , one gets the maximum beam-beam tune shift with respect to 
the order of nonlinear resonance, m, as shown in Fig. |3|, where each maxi- 
mum beam-beam tune corresponds to each dominating multipole resonance. 
For flat beams, it is obvious that if the horizontal tune is not well chosen, 
the Cy,max can De 0.032 instead of 0.0447, however, if the vertical resonances 
have been successfully avoided before £ x< fb reaches its limit, one could possi- 
bly obtain £, x , ma x,FB(jbb = 1 hour)= 0.0632 even difficultly. What should be 
stressed is that in choosing the working point in the tune diagram, one has to 
pay attention to the nonlinear resonances of order as high as 14. To explain 
qualitatively why the maximum beam-beam tune shifts for both round and 
flat beams seem to be limited by the lowest order of resonance, i.e., the 1/4 
resonance, we have plotted in Fig. f| the sum of the multipole strengths from 
m=4 to m=14 assuming that they have the same strength, as expressed: 



m=14 



A(Q y ) = (-l) m/2 sin(27rmQ,) (49) 



m=4 



On the same figure we have plotted also the first term (octupole) in this 
summation with two opposite phases as compare references, and it is ob- 
vious that except two regions of Q y , (0.2 to 0.3) and (0.7 to 0.8), one has 
always that the amplitude of the sum is almost the same as that of the oc- 
tupole term, and in this case the dangerous Q y values are 0.225, 0.275, 0.725 
and 0.775. Another reason for the lowest resonance dominating is that the 
lower the resonance order the more stable the resonance facing to the phase 
perturbations. 

Now we discuss briefly the choice of tunes (working point). Limited to one 
IP and the flat beam case, based on the original work of Bassetti (LNF-135, 



Frascati, Italy), B. Richter has shown in ref. |Z| that the tune should be 



chosen just above an integer or half integer to make a best use of dynamic beta 
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effect, and this conclusion has been experimentally observed in CESR [^3 
Combining this information with what suggested by eq. ^9], one concludes 
that the tune Q should be chosen in the regions (0,0.2) or (0.5,0.7) to obtain 
a maximum luminosity. If the collision is effectuated with a definite crossing 
angle some important synchrobetatron nonlinear resonances, such as 3Q X ± 
Q s = p, should be avoided also. More discussions on the crossing angle 
effects will be given in section 8. Taking CESR and PEP-II for examples, for 
CESR one finds Q x = 10.523 and Q y = 9.597 @, and for PEP-II the actual 
operation working points are Q x = 0.610 ± 0.01 and Q y = 0.580 ± 0.01 for 
Low Energy Ring (LER) and Q x = 0.565 ± 0.01 and Q y = 0.585 ± 0.01 for 
High Energy Ring (HER) |24 , which in principal consist with our suggestion. 

In this paper, under the assumption that the two colliding beams always 
have the same transverse dimensions, we have arrived at the beam-beam 
effect determined lifetimes expressed in eqs. [TJ, and J6". For a given min- 
imum normalized (with respect to the damping time) beam lifetime one gets 
universal maximum beam-beam tune shift values corresponding to different 
cases. In a real machine the situation can be more complicated, such as the 
flip-flop phenomenon which breaks the symmetry assumed above, and in this 



case one can continue the discussion starting from eqs. 25, Eq, and U7\ by 



differentiating a y from a*^, and by replacing a, a x by E/ \/2, S x ./ \/2, respec- 



tively, where £ = yof + cr 2 and S x = J^ + tTj. We will not, however, 
continue our discussions in this direction in an exhaustive way. 



7 Applications to some machines 



Let us look at three machines, PEP-II B-Factory [25| and DA$NE p6[, and 



BTCF |^]], and the first two have been put to operation. The relevant 
machine parameters are shown in Table 2. Figs. [5] and [|give the theoretical 
estimations for the beam-beam limited beam lifetimes in both PEP-II LER 
and HER. Figs. [7] and § show the beam lifetimes versus the beam-beam 
tune shifts in both LER and HER. It is obvious that the nominal charge in 
the bunch of HER is close to the limit which sets the beam lifetime in low 
energy ring, however, the beam lifetime in HER is much longer than that in 
LER. The theoretical results consist with the experimental observation [J25|J . 
Fig. |^ shows the beam lifetime prediction for the DA$NE e + e~ collider with 
single IP. Finally, we study the beam-beam limited beam lifetime in BTCF 
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Machine 


N e 


PyJP cm 


a x j P jum 


dyjp jum 


7 


r y ms 


PEP-II LER 


6xl0 lu 


1.5 


157 


4.78 


6067 


30 


PEP-II HER 


2.8xl0 lu 


1.5 


157 


4.78 


17613 


18.3 


DA$NE 


8.9xl0 lu 


4.5 


2100 


21 


998 


35.7 


BTCF 


1.4xlO u 


1 


450 


9 


3914 


31 



Table 2: Machine parameters 

(standard scheme) and the theoretical result is given in Fig. [K] where the 
dot indicates the designed beam lifetime. 



8 Discussion on the collision with a crossing 
angle 



To get a higher luminosity one could run a circular collider in the multibunch 
operation mode with a definite collision crossing angle. Different from the 
head-on collision discussed above, the transverse kick received by a test par- 
ticle due to the space charge field of the counter rotating bunch will depend 
on its longitudinal position with respect to the center of the bunch which the 
test particle belongs to. In this section we consider first a flat beam colliding 
with another flat beam with a half crossing angle of 4> in the horizontal plane. 
Due to the crossing angle the two curvilinear coordinates of the two collid- 
ing beams at the interaction point will be no longer coincide. The detailed 
discussion about the coordinates transformation can be found in ref. 
When the crossing angle is not too large one has: 



28 



x 



X + zc 



(50) 



where x* is the horizontal displacement of the test particle to the center of the 
colliding bunch, z and x are the longitudinal and horizontal displacements of 
the test particle from the center of the bunch to which it belongs. Now we 
recall eq. [17] which describes the Hamiltonian of the horizontal motion of a 



test particle in the head-on collision mode, and by inserting eq. |50| into eq. 
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|T7| we get: 

1 oo 

(x + z<f)) 8 + ■ ■ ■) Yl S(s-kL) (FB) (51) 



3360a,8 



Since the test particle can occupy a definite z within the bunch according to 
a certain probability distribution, say Gaussian, it is reasonable to replace z 
in eq. ^T] by a z , and in this way we reduce a two dimensional Hamiltonian 
expressed in eq. |5l| into a one dimensional one. What should be noted is that 
eq. |51] takes only the test particle's longitudinal position into consideration 
which is regarded as a small perturbation to the head-on collision case, and 
the geometrical effect will included later. To simplify our analysis we consider 
only the lowest synchrobetatron nonlinear resonance, i.e., 3Q X ± Q s = p 
(where Q s is the synchrotron oscillation tune, and p is an integer) which turns 
out to be the most dangerous one [^J |5D| . Following the same procedure in 
section 4 one gets the dynamic aperture due to the lowest synchrobetatron 
nonlinear resonance as follows: 

A syn . beta , x (s) - [ msip)3 ) NeTeaA ^ 
and 2 

<tj _ A S yn-beta,x( s ) 2 ( _j_ | ( ^\ 

I^syn-beta,, ~ ^ 2 - ^ ^ J {™) 

where $ = —<b is Piwinski angle. Now we are facing a problem of how to 
combine the two effects: the principal vertical beam-beam effect and the 
horizontal crossing angle induced perturbation. To solve this problem we 
assume that the total beam lifetime due to the vertical and the horizontal 
crossing angle beam-beam effects can be expressed as: 

r% + r; ( 1 

T bb,total _ 4 1 1 _| 1 I 

yn — beta,x 

ex P ( i i 1 i — ) ( FB ) ( 54 ) 
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where TZ V} s,fb corresponds to eq. [31|. After the necessary preparations, we 
can try to answer two frequently asked questions. Firstly, for a machine 
working at the head-on collision beam-beam limit, how the beam lifetime 
depends on the crossing angle? Secondly, for a definite crossing angle, to keep 
the beam lifetime the same as that of the head-on collision at the beam-beam 
limit, how much one has to operate the machine below the designed head-on 
peak luminosity? To answer the first question we define a lifetime reduction 
factor: 

= T -^2L (FB) (55) 

T bb,y 



where r£ by is given in eq. ^3], and R(&) will tell us to what extent one 
can increase $. Concerning the second question, one can imagine to reduce 
the luminosity at beam-beam limit by a factor of /($) in order to against 
the additional lifetime reduction term due to the definite crossing angle. 



Physically, from eq. |54] one requires: 



A syn -beta,x(sY\ * + ( A,yna,crossin 9 ,8,y( 




<T v (sf 

dyna,head—on,8,yi,^ \ 



a y (s) 2 



(FB) (56) 



Mathematically, one has to solve the following equation to find the peak 
luminosity reduction factor /($): 



,max,FB 



(FB) (57) 



— bo + Jbn + 4a c 

/(*) = V — (FB) (58) 



where a = 3n 2 £l 4esign FB § 2 /2, b = c = V27c^ maXtFB /3, and £ x , ma x,FB 
0.0447. In fact, /($) corresponds to the luminosity reduction due to the 
synchrobetatron resonance, and to find out the total luminosity reduction 



factor, one has to include the geometrical effects 31j||p2|| . The total luminos- 
ity reduction factor can be expressed as follows: 

F($) = /($) (l + $ 2 )~ V2 (FB) (59) 
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where hourglass effect is no taken into account (i.e. (3 y ,ip > &z)- Tak- 
ing KEKB factory as an interesting example [p5fl , one has a x = 90 /im, 
a z = 0.4 cm, <fi = 11 mrad, $ = 0.49, £, x ,design ) FB = 0.039, and by putting 
$ = 0.49 into eq. [58] one finds F(0.49) = 83.5% which is very close to a 
three dimensional numerical simulation result, i.e., 85% of designed head-on 
luminosity, given in ref. [33|. In Figs. [11] and [12] we show how -R($) and 
F($) depend on Piwinski angle where £, x ,design,FB = 0.039 has been used in 
Fig. [12|. Finally, when the crossing angle is in the vertical plane or the beam 
is round, one gets: 



K 



syn—beta,y 



3tt 2 v e$ 



(FB) 



(60) 



and 



where r 
|54] by eq. |60 



K 



syn—beta,y 



32 



1 



277T 2 i&$ 



(RB) 



(61) 



(Ty/(T X and $ = ^0 as defined before. Replacing TZ syn -beta,x in eq. 

and following the same procedure shown above one can 



or eq. 



easily make the corresponding discussions on the luminosity reduction effects. 
What should be remembered is that the geometrical luminosity reduction 
factors for the vertical crossing angle and the round beam cases are (1 + 
($/ r ) 2 )-V2 anc } _j_ $2^ respectively. 



9 Conclusion 

In this paper we have established analytical formulae for the beam-beam 
interaction limited dynamic apertures and beam lifetimes in e + e~ circular 
colliders for both round and flat beam cases. It is shown analytically why for 
flat colliding beams one has always £, y , max around 0.045 and why this value 
can be almost doubled by using round beams. Applications to the machines, 
such as PEP-II, DA$NE, and BTCF have been made. Finally, the luminos- 
ity reduction due to a crossing angle has been discussed, and an analytical 
formula for the luminosity reduction factor is derived and compared with a 
numerical simulation result for KEKB factory. 



16 



10 Acknowledgement 



The author thanks J. Haissinski for his careful reading of the manuscript, 
critical comments, and drawing my attention to the paper of B. Richter. 
Stimulating discussions with A. Tkatchenko and J. Le Duff are appreciated. 



References 

[1] J. Gao, "Analytical estimation of the dynamic apertures of circular ac- 
celerators", Nucl. Instr. and Methods, A451 (3) (2000), p. 545. 

[2] A. Piwinski, "Observation of beam-beam effects in PETRA", IEEE 
Trans, on Nucl. Science, Vol. NS-26, No. 3, June 1979. 

[3] M.A. Furman, "Beam-beam issues for high luminosity e + e~ colliders", 
Frascati physics series Vol. X (1998), p. 123. 

[4] E. Keil, "Beam-beam dynamics", CERN 95-06, p. 539. 

[5] R. Talman, "Resonances in accelerators", AIP Conference Proceedings 
153 (1984), p. 835. 

[6] T. Chen, J. Irwin, and R.H. Siemann, "High-order horizontal resonances 
in the beam-beam interaction", Nucl. Instr. and Methods, A402 (1998), 
p. 21. 

[7] S. G. Peggs and R. M. Talman, "Nonlinear problems in accelerator 
physics", Ann. Rev. Nucl. Part. Sci. Vol. 36 (1986), p. 287. 

[8] S. Peggs and R. Talman, "Beam-beam luminosity in electron-positron 
colling rings", Physical Review D, Vol 24, No.9 (1981), p. 2379. 

[9] A. Chao and M. Tigner (Editors), "Handbook of accelerator physics and 
engineering", World Scientific (1999), p. 134. 

[10] M. Bassetti and G. Erskine, "Closed expression for the electrical field of 
a two-dimensional Gaussian charges", CERN/ISR-TH/80-06 (1980). 

[11] V. Ziemann, "Beyond Bassetti and Erskine: beam-beam deflections for 
non-Gaussian beams", SLAC-PUB-5582 (1991). 



17 



[12] R. Talman, "Multiparticle phenomena and Landau damping", AIP Con- 
ference Proceedings 153 (1984), p. 790. 

[13] K. Hirata, "The beam-beam interaction: Coherent effects", AIP Con- 
ference Proceedings 214 (1990), p. 175. 

[14] K. Hirata, "Coherent betatron oscillation modes due to beam-beam in- 
teraction", Nucl. Instr. and Methods, A269 (1988), p. 7. 

[15] K. Hirata and E. Keil, "Barycentre motion of beams due to beam- 
beam interaction in asymmetric ring colliders", Nucl. Instr. and Meth- 
ods, A292 (1990), p. 156. 

[16] S. Krishnagopal and R. Siemann, "A comparison of fat beams with 
round", AIP Conference Proceedings 214 (1990), p. 278. 

[17] S. Krishnagopal and R. Siemann, "Simulation of round beams", PAC89, 
p. 836. 

[18] M. Bassetti and M.E. Biagini, "A beam-beam tune shift semi-empirical 
fit", Frascati physics series Vol. X (1998), p. 289. 

[19] J. Gao, "Analytical formulae for the maximum beam-beam tune shift in 
electron storage rings", Nucl. Instr. and Methods, A413 (1998), p. 431. 

[20] J. Gao, "Analytical formulae for the maximum beam-beam tune shift 
in circular colliders", LAL-SERA-99-148. 

[21] H. Bruck, "Accelerateurs circulaires de particules" (Press Univer sit aires 
de France, 1966, p. 268. 

[22] B. Richter, "Design considerations for high energy electron-positron 
storage rings" , Proceedings of the international symposium on electron- 
positron circular colliders", 26-30 September, 1966, Saclay, France, p. 
1-1-1. 

[23] D. Sagan, "The dynamic beta effect in CESR", Proceedings of PAC95, 
USA (1995), p. 2889. 



18 



[24] M. Placidi, "Beam-beam issues from the recent PEP-II commissioning", 
Proceedings of workshop on beam-beam effects in large hadron colliders, 
CERN, 1999, p. 45. 

[25] J.T. Seeman, "Commissioning results of the KEKB and PEP-II B- 
factories", PAC99, New York, 1999, p. I. 

[26] G. Vignola, "DA^NE, the Frascati ^-factory" , PAC93, USA (1993), 
p. 1993. 

[27] "Feasibility study report on Bejing Tau-Charm Factory", IHEP-BTCF 
Report-01, Dec. 1995. 

[28] K. Hirata, "Analysis of beam-beam interactions with a large crossing 
angle", Physical Review Letters, Vol. 74, No. 12, March 1995, p. 2228. 

[29] A. Piwinski, "Satellite resonances due to beam-beam interaction" , IEEE 
Trans, on Nucl. Set., Vol. NS-24, No. 3, June 1977, p. 1408. 

[30] A. Piwinski, Computer simulation of the beam-beam interaction at a 
crossing angle", IEEE Trans, on Nucl. Sci, Vol. NS-32, No. 5, Oct. 
1985, p. 2240. 

[31] A. Piwinski, "Satellite resonances due to beam-beam interaction", Nucl. 
Instr. and Methods, 81 (1970), p. 199. 

[32] N. Toge and K. Hirata, "Study on beam-beam interactions for KEKB" , 
KEK Preprint 94-160. 

[33] K. Ohmi, "Simulation of the beam-beam effect in KEKB", KEK pro- 
ceedings 99-24, p. 187 (International workshop on performance improve- 
ment of electron-positron collider particle factories). 



19 



1000000 



100000 r 



10000 



1000 



100 



10 



0.04 0.05 0.06 0.07 0.08 0.09 0.1 
Vertical beam-beam tune shift 



T 

Figure 1: vs £* for flat beam case. 
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Figure 2: ^ vs £* for round beam case. 
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Figure 3: The maximum beam-beam tune shift vs the resonance order. For a 
flat beam the dotted line and the diamonded line correspond to the horizontal 
and vertical tune resonances, respectively. The triangled line corresponds to 
the vertical resonance for a round beam case. 
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Figure 4: The nonlinear perturbation amplitude vs the vertical tune. The 
two inversely phased sinusoidal curves correspond to the amplitude of the 
octupole term, and the fast oscillating curve is the sum of the multipoles of 
the same amplitude from octupole to 28 poles. 
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PEP-II Low Energy Ring 




Figure 5: The beam lifetime due to the beam-beam interaction vs the particle 
population in the bunch in the low energy ring of PEP-II. 
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Figure 6: The beam lifetime due to the beam-beam interaction vs the particle 
population in the bunch in the high energy ring of PEP-II. 
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Figure 7: PEP-II Low Energy Ring (flat beam): the beam-beam limited 
lifetime r 6 * 6 (hours) vs the vertical beam-beam tune shift £ y . 
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Figure 8: PEP-II Low Energy Ring (round beam is assumed): the beam- 
beam limited lifetime r b * fe (hours) vs the vertical beam-beam tune shift £ r 
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Figure 9: The beam lifetime due to the beam-beam interaction vs the particle 
population in the bunch in the DA$NE ^-Factory. 
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Figure 10: The beam lifetime due to the beam-beam interaction vs the par- 
ticle population in the bunch in the Beijing r-C Factory (standard scheme). 
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Figure 11: The lifetime reduction factor vs Piwinski angle 
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Figure 12: The luminosity reduction factor vs Piwinski angle $. The 

curve is obtained by taking £, X) design,FB = 0.039 (KEKB), and the dot is the 
numerical simulation result given in ref. [|33 
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